Upward curvature of the upper critical field 
in the Boson— Fermion model 
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We report on a non-conventional temperature behavior of the upper critical field [HciiT)) which 
is found for the Boson-Fermion (BF) model. We show that the BF model properly reproduces two 
crucial features of the experimental data obtained for high-Tc superconductors: IIc2{T) does not 
saturate at low temperatures and has an upward curvature. Moreover, the calculated upper critical 
field fits very well the experimental results. This agreement holds also for overdoped compounds, 
where a purely bosonic approach is not applicable. 



Some of the unusual properties observed experimen- 
tally in the high temperature superconductors (HTSC) 
can be explained in terms of the effective two component 
models, which involve boson and fermion degrees of free- 
dom. There are numerous examples of such theoretical 
scenarios, for instance: electron fields coupled to gauge 
fluctuations, the RVB spinon-holon theory, the coupled 
electron-phonon systems, etc. In this paper we consider 
one of these models, where the conduction band parti- 
cles (fermions) coexist and interact with the localized 
electron pairs (hard-core bosons). This, so called Boson 
Fermion (BF) model, had initially been introduced a 
couple of years before the discovery of HTSC materials. 
The authors have conjectured that it might describe an 
emerging physics of the electron-phonon systems for the 
intermediate coupling strength. 

From a formal point of view, the Hamiltonian of this 
BF model can be considered as an effective one for several 
models, which are often used in the solid state physics. It 
has been derived so far: (1) from the generalized periodic 
Anderson model with the large on-site attraction by elim- 
inating the hybridization between the wide and narrow 
band electrons via canonical transformation [||, (2) from 
a purely fermionic single band system described by the 
extended Hubbard model using the concept of bosoniza- 
tion for lattice fermions (3) from the Hubbard model 
on the two dimensional plaquettized lattice in the strong 
interaction limit using the contractor method Q| , and (4) 
from the resonating valence bond state of the t — J model 
in the path integral technique 

The same BF model has been considered also by many 
other authors who postulated it in a more ad hoc way. 
For example in the Ref. Q author has proposed the BF 
type charge transfer Hamiltonian basing his arguments 
on interpretation of the optical experiments. Other ideas 
could be found in the Ref. . Authors have represented 
the patches of the 2D Brillouin zone near the so called 
hot spots via dispersionless bosons, which are coupled 
to fermions of a remaining part of the Brillouin zone. 



The BF model is not only considered in a context of 
HTSC. There are recent attempts to apply the same type 
of picture for a description of the magnetically trapped 
atoms of alkali metals |^ . 

Unconventional mechanism of superconductivity 
within the BF model has been studied in a number of 
papers |l|, |, 0, [lO[ |ll|, |l|, |l|, [l4|. Fermions acquire 



superconducting coherence via exchange of the hard-core 
bosons. The same processes lead also to the effective 
mobility of bosons. Fermions/bosons can undergo a 
phase transition into the superconducting/superfluid 
phase at identically the same critical temperature 
Tc = Tf^ = T^E @- It is worth to point out some 
of the unusual properties obtained for the BF model 
which are known in the most of the HTSC materials: 
{i) a non-BCS ratio As{T = 0)/kTc > 4 (except for 
the far under- and over-doping regimes) |l^; (m) 
linear in T resistivity of a normal phase up to very 
high temperatures ||]; (iii) change of sign of the Hall 
constant above Tc and the anomalous Seebeck coefficient 
0; (iv) appearance of the pseudogap in a normal 
phase for temperatures T* > T > (v) a 

particle-hole asymmetry of the single particle excitation 
spectrum in the normal phase |9| . 

Some of unusual properties of HTSC are related to 
the upper critical field. Hc2 can achieve values of a 
few hundred Tesla. Moreover, the resistivity measure- 
ments for HTSC clearly show an upward curvature of 
Hc2(T), with no evidence of saturation even at low tem- 
peratures From the theoretical point of view 
the upward curvature of the critical field occurs for in- 
stance in: Bose-Einstein condensation of charged bosons 
p2| , Josephson tunneling between superconducting clus- 
ters and in mean-field-type theory of Hc2 with a 
strong spin-flip scattering However, this feature 
cannot be explained within a conventional theory of Hc2 
Ipsf . Therefore, it is a natural test for theoretical ap- 
proaches to HTSC. 

In the present study we show that the upward curva- 
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ture of Hc2{T) is an intimate feature of the BF model. 
This result strongly supports the BF Hamiltonian as a 
model of HTSC. 

We consider the Hamiltonian of the two-dimensional 
BF system, immersed in a perpendicular, uniform mag- 
netic field 



+ 5Z - blh + vY, {bUic,l + h.c.) (1) 

i i 

We use the standard notation for annihilation (creation) 
operators of fermion Cia (cj^.) with spin a and of the hard 
core boson bi (6|) at site i. Fermions interact with bosons 
via the charge exchange interaction v. fj, is the chemical 
potential and tij (A) is the hopping integral that depends 
on the magnetic field through the vector potential A: 

(A) = (0) ^ A • dl^ 

To proceed, we apply the mean field decoupling for the 
boson fermion interaction 



(2) 



which is justified when v is small enough in comparison 
to the kinetic energy of fermions. After the decoupling 
we deal with the effective Hamiltonian composed of 
the fermion and boson contributions H ~ + H^, 
coupled through the selfconsistently determined expec- 
tational values (cijCif ) and (bi) |l^ 

^ [ty(A) - S,jn] cl^Cj^ + (p*c,^c,t + h.c.) , 



H- ^Y 



[iAB~2^l) blb, + A, bl + A*b, 



(3) 
(4) 



where Ai = v (ciiCii) and pi = v {bi). 

One can exactly diagonalize the bosonic subsystem us- 
ing a suitable unitary transformation. Statistical expec- 
tation values of the number operator blbi and the param- 
eter Pi are given by |0, O 



blb^ = i-^^-^tanhm. 



(5) 
(6) 



where = \ \/{Ab - 2p)^ + 4|Ai|2 and k is the Boltz- 
mann constant. At the phase transition the supercon- 
ducting order parameter is infinitesimally small and one 
can expand pi in powers of (ciiCi^) up to the leading 



order. Then, the fermionic subsystem is described by 

- V(.T)Y,((cl^cl)c.ic,^+h.c), (7) 



where 



ViT) 



Ab-2h 



■ tanh 



A, 



2p 



2kT 



(8) 



The effective Hamiltonian (0) represents a BCS-type 
model with isotropic on-site pairing. The boson-fermion 
coupling enters through the pairing interaction (|^) 
and the chemical potential. The latter quantity should be 
evaluated from the conservation of the total charge ntot = 
2ns +nF ^ 2/NJ:^ ib\h) + (4^..). We 

restrict further investigation only the nearest-neighbor 
hoping, when the fermionic energy spectrum is known as 
the Hofstadter butterfly H). 

First, we assume that the bosonic level is in the middle 
of the fermionic band {As =0) and ntot = 2. It is the 
simplest case, when the parameters of the Hamiltonian 
can be evaluated without a numerical investigation 
of the Hofstadter butterfly. Namely, since the Hofstadter 
energy spectrum is symmetric, for /i = one gets n^ = 1, 
ub = 1/2, ntot = 2 and V(T) = v^/{2kT). In order to 
calculate Hc2 one can apply a two-dimensional version 
of the Helfand-Werthamer (HW) theory where the 
coupling constant, A, depends on temperature. Since, 
A — V p-ps, where p-ps is the density of states at the Fermi 
level, one gets 



A(T) = A(r,) 



T' 



(9) 



The HW theory was derived for a free-electron gas and 
neglects the Landau level structure (so called quasiclas- 
sical limit). However, for a weak magnetic field these 
effects do not lead to an essential modification of Hc2 (T) 
Figure 1 shows Hc2(T) calculated within the HW 
approach with the coupling constant determined by Eq. 
(^. In order to predict the qualitative temperature de- 
pendence of the critical field it is sufficient to recall the 
following facts: (i) when A(r) = const (standard HW 
theory), Hc2{T) is linear for a weak magnetic field; (ii) 
in the BF model A increases with the decrease of tempera- 
ture and diverges when T — > (Eg. |^). Then, it becomes 
obvious that BF model properly reproduces two crucial 
features of the experimental data: Hc2{T) does not sat- 
urate at low temperatures and has an upward curvature, 
at least for a weak magnetic field. Numerical results (see 
Fig. 1) show that the latter property holds in the whole 
range of temperature. 

Next, we show that the BF model properly describes 
Hc2{T) for a wide range of model parameters, when the 
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FIG. 1: Temperature dependence of the upper critical field, 
obtained from the IfW approach with the coupling constant 
given by EqM . We have assumed kTc{H = 0) = 0.02t. 



problem cannot be reduced to the effective HW theory. 
In this case we apply a lattice version of the Gor'kov 
equations 



-UJn 



(10) 



Here, G{i^j,u)n) is the one-electron Green's function in 
the presence of a uniform and static magnetic field and 
UJn is the fermionic Matsubara frequency. With the help 
of the Hofstadter approach |2^, equation (|l^) can ex- 
actly be solved for clusters of the order of lO"' lattice 
sites. For the details we refer to Ref. [|2^. In contradis- 
tinction to the quasiclassical approaches (e.g., HW or 
Ginsburg-Landau theory) we explicitly account for the 
actual structure of Landau levels. 

Figure 2 shows the upper critical field obtained for 
ntot = 1 and different positions of the bosonic level. We 
consider two cases: (z) when this level is below the Fermi 
energy (As < 0) there is a finite number of bosons also at 
T —f (m) for As > bosonic states are occupied only 
virtually. The upward curvature of Hc2iT) appears pre- 
dominantly in the first case (solid line in Fig. 2). When 
As is shifted above the Fermi energy, the curvature is 
gradually reduced. Finally, when A^ 3> kT the curva- 
ture changes from positive to negative (dashed line in 
Fig. 2) and one reproduces standard results for a purely 
fermionic system |2^ . In Fig. 2 we have not presented re- 
sults in the low temperature regime. Although, we have 
carried out numerical calculations for large clusters, this 
approach is not applicable at genuinely low temperatures. 
In this case the Cooper pair susceptibility accounts only 
for very few fermionic states (with energy close to the 
Fermi level) instead of a continuous density of states. 
However, one can prove that BF model qualitatively re- 
produces Hc2{T) also for T ^ 0. Combining Eqs. (||) 
and (P) one can express the effective pairing potential in 



FIG. 2: Hc2(T) obtained from Eq. (jig) for ntot = 1 and 
different bosonic levels. The fermion-boson coupling was ad- 
justed to obtain the same critical temperature in the absence 
of magnetic field. 




FIG. 3; Fit to experimental results (continuous curve). The 
experimental data (points) were taken for Tl2Ba2Cu06 from 
Ref. [pi, and for Bi2Sr2CuOy from Ref. H. 



terms of the bosonic occupation number 
2i;2 



V{T) 



2fi 



bib, 



(11) 



It is straightforward to note that for {h\hij — 1/2 one 
gets As = 2^ and Eq. (|) is reduced to V{T) = v^jlkT. 
On the other hand, when lb\bi^ ^ 1/2 but < ub < 1 
the denominator in Eq. ( pT| ) vanishes for T ^ 0. There- 
fore, V{T) diverges for T ^ provided that < Ub < 1. 
It means that a requirement of partial occupation of 
bosonic states is sufficient to reproduce the experimental 
low temperature behavior of i?c2 • 

To complete the discussion, we show that BF model 
accurately reproduces the experimental data (see Fig. 
3). We have chosen an appropriate set of model param- 
eters, for which Hc2{T) fits very well the results pre- 
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sented in Refs. [EOl 21 . The theoretical curve was cal- 



culated for bosonic level that is situated slightly above 
the Fermi energy (ntot = 1 and =0). In this case 
Up ^ ub and the BF system is mainly of fermionic char- 
acter. The experimental data were obtained for over- 
doped compounds, when HTSC exhibits a Fermi liquid 
character Therefore, BF model is fully consistent 

with fermionic-type behavior of the overdoped HTSC, 
whereas the purely Bosonic approach is not applica- 
ble in this regime. 

To conclude, we have shown that the BF model ac- 
curately describes the upper critical field observed in 
HTSC. It is known that the BF model correctly describes 
the pseudogap phenomenon and some other unconven- 
tional normal-state properties of HTSC. In this letter 
we have shown that this model reproduces also other un- 
usual feature of HTSC, which cannot be explained within 
a standard BCS-type approach. Moreover, our results 
remain in agreement with the Fermi liquid behavior of 
the overdoped compounds. Therefore, it is an important 
support for the BF Hamiltonian as a model of HTSC. 
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